Vortices in superfluid trapped Fermi gases at zero temperature 
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We discuss various aspects of the vortex state of a dilute superfluid atomic Fermi gas at T = 0. 
The energy of the vortex in a trapped gas is calculated and we provide an expression for the 
thermodynamic critical rotation frequency of the trap for its formation. Furthermore, we propose 
a method to detect the presence of a vortex by calculating the effect of its associated velocity field 
on the collective mode spectrum of the gas. 



Inspired by the impressive progress in recent years in 
the field of Bose-Einstein condensation in dilute atomic 
gases, increasing attention is being devoted to examine 
the behavior of a gas of fermionic atoms at the same 
ultra- low temperatures. Experimentally, the trapping 
and cooling of fermionic alkalis has been demonstrated 
reaching temperatures as low as ~ Tp/4 for '"'K and 
^Li with Tp denoting the Fermi temperature. It is 
well known from condensed matter and nuclear physics 
that a gas composed of two different internal states of the 
same fermionic particle which interact via an attractive 
interaction is unstable to formation of so-called Cooper 
pairs, thus becoming a superfluid. After the possibility 
of such superfluid transition for trapped Fermi gases was 
proposed ||^, a lot of theoretical work has been focus- 
ing on various properties of this system ^ . At the same 
time a major experimental goal has become to observe 
the formation of the superfluid state. 

One of the intriguing properties of a superfluid is the 
possibility of forming quantized vortices. For a Bose- 
Einstein condensate, the study of vortices has produced 
several interesting results 0. Recently, some aspects of 
the vortex state of a trapped superfluid Fermi gas close to 
the critical temperature Tc of the superfluid phase tran- 
sition were considered Q . In this paper we are interested 
in the properties of the vortex state of clouds of trapped 
Cooper-paired fermions at T = 0, and in particular in un- 
derstanding under which conditions a vortex forms, what 
is its energy, and how it can be detected. We consider 
large systems where the coherence length ^ of the super- 
fluid is much smaller than the extent of the cloud. In this 
limit, we are able to derive an analytical estimate of the 
energy of a vortex in a trap, thereby predicting the criti- 
cal rotation frequency for its formation. Also, we propose 
a way of observing the vortex by calculating its effect on 
the collective mode spectrum of the gas. The paper is 
organized as follows: flrst in Sec. ||, we examine for which 
values of the characteristic parameters the vortex is well 
localized within the gas. In Sec. || we present a simple 
model for calculating the energy of a vortex in a uniform 
superfluid Fermi system. Using the result of Sec. |l|, in 
Sec. [II we calculate the energy of a singly quantized vor- 
tex in a trapped gas and from that we obtain the value 
of the thermodynamic critical rotation frequency for its 
formation. The problem of observing the vortex state is 



considered in Sec. IV, where we calculate how the pres- 



ence of a vortex influences the collective mode spectrum 
of the gas. Finally, we summarize our results in Sec. 
Given the uncertainties intrinsic in any simple model of 
the vortex, such as the one presented in Sec. |l|, in App. 

we briefly discuss another possible way of describing a 
vortex in a uniform gas, based on an approximate zero- 
temperature Ginzburg-Landau approach We then 
calculate what this alternative method gives for the en- 
ergy of the vortex, and compare the two results to show 
that they do not differ in any significant way. 



I. BASIC CONSIDERATIONS 

In the dilute ultracold limit the effective interaction 
between identical fermionic atoms vanishes due to the 
Pauli principle, and that between different ones can be 
well described by one parameter only, the s-wave scat- 
tering length a. For a negative scattering length, the 
interaction is attractive and if the number of particles 
in the two internal states is the same the T = ground 
state of the gas is a superfluid. The critical temperature 
Tc for the transition to such a superfluid state in a dilute 
gas was first determined for a uniform system by Gorkov 
and Melik-Barkhudarov |p^ , and using a more modern 
approach by Heiselberg et al. The predicted value 

is 
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where A stands for 2kp\a\/iT, ep is the Fermi energy com- 
mon to the two species of fermions, kp the associated 
Fermi wavenumber, and 7 ~ 1.781 is related to Euler's 
constant C by 7 = e*-^. The pairing gap A at T = is, as 
usual in BCS theory, related to the critical temperature 
by Ao=TTj-^kBTc lllJl|. 

When applying this result to a gas trapped by a har- 
monic oscillator potential, as in the cases of experimental 
interest today, some requirements have to be met. The 
first one is, just as for the uniform case, that the density is 
everywhere so low that the gas is dilute, i.e. fci?(r)|a| <C 1. 
We have introduced a local Fermi wavenumber kp(r). 
This corresponds to using the Thomas-Fermi approxi- 
mation, which is valid if ep ^ huix, where ujt is the 
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frequency of the oscillator (which for the time being we 
assume to be isotropic). This condition is always satis- 
fied if the particle number is sufficiently large, since for 
a harmonic potential ep = {GNcrY^^hujT, with TVo- being 
the number of particles of one species. 

Another condition for applicability of Eq. (0) is that 
ksTc ^ hujT 1^5. When this latter condition is not 
satisfied, the shell structure of the harmonic oscillator is 
crucial when determining the superfluid properties of the 
gas, and Eq. ([|) in general breaks down. 

In a superfluid Fermi gas at zero temperature the size 
of the vortex core of a singly quantized vortex is given 
approximately by the BCS coherence length ^bcs = 
hvp /ttAq, where vp = Tikp/ma is the Fermi velocity 
and nia the mass of a single atom. It is clear that 
in order for a vortex to appear at all, the BCS coher- 
ence length (size of the vortex core) at the center of the 
cloud has to be smaller than the size of the cloud itself, 
which in the Thomas-Fermi approximation is given by 
Rtf = {2ep/maUj'^Y^'^- If this were not so the super- 
fluid properties of the system would be more like those 
of a nucleus (for which ^ > i?) than those of a bulk su- 
perfluid. 

Substituting the appropriate expressions one can im- 
mediately see that ^bcs/Rtf = n'^hLUT/^Q- So that 
requiring S,bcs ^ Rtf corresponds to demanding that 
Ao 3> TT^^hujT- This condition is automatically satisfied 
if ksTc 3> hujT, but is not at all obviously realized in 
possible practical circumstances. Indeed, if we assume 
the validity of equation (0), and of the related value of 
Aq, and we use the expression ep = {GN^Y^^hujT for the 
Fermi energy, we obtain 
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We may then immediately see that unless A, and there- 
fore fci?|a|, is sufhciently close to one, the exponential is 
very large and the condition in Eq. (^ is not satisfied, 
implying that the coherence length is much larger than 
the radius of the cloud and the rotation pattern very dif- 
ferent from a vortex state. If however kp\a\ is too close 
to one (~ 0.3 — 0.4 or more) the formula becomes unre- 
liable because the gas is no more dilute and effects due 
to induced interactions, which strongly modify the value 
of Aq obtained in the dilute limit, must be taken into 
account For the sake of the present work we shall 

not consider these effects; we study regions of densities in 
which Eq. (|l]) is reasonably reliable keeping kp\a\ < 0.4. 
There is then a region of applicability of Eq. (|l|) for a 
trapped gas which depends on the number of particles 
and the scattering length. In order to find this region we 
impose the equality in Eq. (H) , and we plot in figure |l| the 
critical number of atoms N„^c for which £,bcs / Rtf = 1, 
as a function of fc_F|a|. Well above the curve we are in the 
regime where the local density approximation can be ap- 
plied and a vortex may form, and below it the superfluid 
has a character more related to that of a nucleus. Since 



the value of kp\a\ can be simply increased by keeping 
the number of particles fixed and tightening the external 
trapping potential, we see that if No- is sufficiently large 
(^ 10^) these systems have the interesting possibility of 
going from one regime to the other. 

In the remainder of this work we assume that we are 
in the upper region of fig. ^ and therefore that ^bcs ^ 
Rtf- In this region a vortex forms in the cloud if it 
is stirred at an angular velocity greater than a critical 
one LUci , which we shall calculate using a thermodynamic 
approach. 



II. VORTEX IN A UNIFORM GAS 

Let us for the time being suppose that the vortex we 
want to describe is in a uniform gas. In particular we may 
take the system to be in a cylinder of radius Rc ^ S.BCS- 

Associated with the vortex there is a superfluid ve- 
locity flow which decreases with the distance from the 
vortex axis: 'Vy{p) — Kh/2maP, where k is the number 
of quanta of circulation of the vortex. In a simple model 
this velocity field extends from p ^ k£_bcs to p — Rc- At 
distances shorter than ~ k^bcs i the kinetic energy asso- 
ciated with the rotation becomes high enough to break 
the Cooper-pairs, and thus the fluid inside a cylinder of 
radius ~ k^bcs about the vortex axis can be thought of 
as being in a normal (non-superfluid) state. The energy 
per unit length associated with a vortex is then given by 
the sum of two contributions. One is the kinetic energy 
due to the flow: 
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and the other one is the loss in condensation energy about 
the vortex axis 
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where econd = SAqU^/Acp is the condensation energy 
per unit volume due to the pairing and the usual 
expression for S,bcs has been employed. Notice that we 
have here introduced the one-species particle density riu, 
and we have supposed that this is a constant throughout 
the system, since contrary to the boson case it is not the 
particle density but onl y t he pairing field that changes 
close to the vortex axis . 

The total energy per unit length of a vortex is therefore 
given in the simple model by 
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The important feature of this result is to point out that 
for large systems (i.e. for which i?c ^ k£,bcs) the most 
relevant contribution is the logarithmic one arising from 
the kinetic integration. The value of the constant inside 
the logarithm will depend on the choice of the model used 
to describe the vortex. A more reliable value would be 
obtained from a numerical solution of the Bogoliubov-de 
Gennes equations, although it is unlikely that it will differ 
significantly from the one found here, since one expects 
it in any case to be of order one. As an example of what 
a different approach may yield, in App. A we show the 
result for the total energy of the vortex obtained using a 
zero-temperature Ginzburg-Landau model. As we shall 
see one obtains, as foreseen, the same expression as in 
Eq. (j^), with coefficient 1.65 instead of 1.36 inside the 
logarithm. 

For what follows we shall not need to know the precise 
value of this coefficient, which may be better determined 
in the future, and we shall therefore leave it unspecified 
and state our result as 
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with the understanding that D is some contant of order 
one. From Eq. (^) it is already clear that one vortex 
with K = K ^ 1 has greater energy than k vortices with 
K = 1 since in any case we need to have k^bcs ^ Rc- 
This implies that vortices with k ^ I are unstable [ pT[ . 
Therefore for the considerations that follow we shall set 
K = 1. 

With this solution, recalling that the thermodynamic 
critical velocity for formation of a first vortex is given 
by Wci = £v/J~-v [@, and using the fact that the total 
angular momentum per unit length of the system with a 
vortex is ~ fiirR^n^, corresponding to Ti per Cooper 
pair, we can immediately state what the critical velocity 
is in a uniform system, which is of course a well known 
result 
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The result should be compared with the critical velocity 
found for a Bose-Einstein condensate. This is completely 
analogous if the mass of a single bosonic atom is replaced 
with that of a Cooper pair (2ma), and the boson coher- 
ence length by the BCS one. 

Notice that since ^bcs Aq ^, from the measurement 
of LOci one could in principle deduce the value of Ao if D is 
known. This possibility is usually lost however in a non- 
uniform system since several values of Aq are integrated 
over. 



III. VORTEX FORMATION IN TRAPPED GASES 

In this section we calculate the energy of the vortex in 
a trapped gas, specializing to the trapping configurations 
used in a typical experiment. 



The atoms are generally confined in a cylindrically 
symmetric harmonic potential of the form 



(8) 



and the density profile of the gas is, within the Thomas- 
Fermi approximation, given by 



\ 9 2 1 2 \ 3/2 



Rl 



(9) 



Here ricfl — ^^.(O, 0) is the density at the center of the 
cloud and we have taken the profiles of the two species 
to be identical. The anisotropy of the trap is controlled 
by the coefficient At. The energy of a cloud with a vor- 
tex along the z axis can be calculated with the proce- 
dure devised by Lundh et al. One can divide the 
cloud in vertical slices of height dz and use the result 
for a cylinder of radius pi such that £,bcs ^ Pi <C 
R± = Rz/Xt, within which one can assume that the gas 
is approximately uniform. The energy per unit length 
associated with the vortex in a slice at z is then given by 
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where R±{z) = (1 — z^/R'^.^^^Rz/Xt is the value of p up 
to which the cloud extends for a given z, and na-{0, z) is 
the density on the z-axis at height z. The second term in 
Eq. ([l^) gives the kinetic energy of the superfluid outside 
the cylinder of radius pi. 

With n^{p, z) given by (||) we then get 
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This result differs from the boson case in ref. in the 
power 3/2 instead of 1 in the density distribution. Using 
the fact that 
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and that unless z is very close to i?z one can assume 
pi <C Rj_{z), we finally obtain 
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In order to proceed with the z integration we need to 
know the expUcit dependence of (,bcs on z. In the dilute 
gas approximation where Eq. (p|) is vahd this is given by 
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where /cf.o — {'^"maep /'h^Y^'^ and Aq — 2fci?,o|a|/7i' are 
the local Fermi wave-number and A respectively, evalu- 
ated at the center of the cloud. Inserting this value into 
Eq. (p!^), using the expression for R^{z) and integrating 
over z we get after some cumbersome but straightforward 
calculations 
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Note that Eq. (|15|) predicts the energy cost of the vortex 
to be negative for small A:i?.o|a| and /fiuj±_ not too large. 
This is clearly an unphysical result reflecting the fact that 
in the limit of relatively few particles trapped and small 
fcF,o|«J, the condition £,bcs ^ Ri- is violated making 
Eq. ( p^ ) invahd. In the regime £,bcs < R^-, Eq. ( p^ ) 
yields positive vortex energies as expected. If we ignore 
the non-rotating particles at the core of the vortex, the 
total angular momentum of the vortex state is Ly = N^Ti 
and the critical rotation frequency Wci = E^j/L^ for the 
formation of a vortex in a trap given by Eq. (0) is 
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with l± = y''h/'muj± being the harmonic oscillator length 
in the radial direction. For realistic parameters of the gas, 
this critical frequency is rather small: choosing for D the 
value obtained in App. taking fc^.olal = 0.4, uj± — 
ujz = LOT (i.e. At — 1) and ep — 200 Ticjt corresponding 
to isotropic trap with N„ ^ 1.3 x 10^, we obtain cj^i — 
0.0035 oj^. The reason for the critical frequency being 
so small is that the angular momentum per atom is h/2 
yielding Ly = N^Ti whereas the energy given by Eq. ( p5| ) 

2/3 

only scales as Na ■ The Fermi pressure expands the 
cloud and reduces the density at the center of the trap. 
Since the energy of the vortex mainly comes from regions 
close to the vortex axis where the superfluid velocity is 
high, the energy of the vortex is correspondingly reduced. 



IV. OBSERVATION OF THE VORTEX 

Contrary to the situation for Bose-Einstein conden- 
sates, the presence of a vortex in the Fermi gas does 
not alter the density profile significantly ||l^. One can- 
not therefore observe the vortex simply by looking at the 
density profile. It has been suggested to use the laser 



probing method of Ref. Ig^l to detect the local decrease 
of the pairing near the center of the vortex . Here we 
examine a different method based on measuring the col- 
lective mode spectrum of the gas. In the case of no vortex 
present, excitations of the gas carrying equal and oppo- 
site angular momentum along the z-axis are degenerate 
in energy. The velocity field associated with a vortex 
aligned with this axis lifts the degeneracy since the rota- 
tional symmetry is removed; the velocity flow of the exci- 
tation is either parallel or anti-parallel to that of the vor- 
tex giving rise to an energy splitting of the modes [pT|-p3[ . 
Since the collective mode frequencies of the gas can be 
measured with a fairly high precision, the possibility of 
detecting the presence of the vortex by its spectroscopic 
signatures is a promising method. Indeed, this method 
has proven to be very useful in the case of a vortex in 
a BEC j2^. The calculations will be carried out for an 
isotropic trap with R^ = Rz — Rtf and uj± — lu^ — ujt- 
In the £,Bcs ^ Rtf limit considered in this paper, the 
collective modes of the superfluid gas for T = can be 
calculated using a hydrodynamic theory. The relevant 
continuity and superfluid velocity equations read p5| : 



dtn{r,t) = -Vj(r,i) 

9tv,(r,i) = —W[mavl/2 + HF 

m„ 



(17) 



with ns(r,t), n„(r,i), and n(r,t) — ns{r,t) -f n„(r,i) 
being the superfluid-, normal- and total density of 
the gas respectively. The total current is j(r, i) — 
Us (r , i) (r , t) -I- n„ (r , i) v„ (r, t) , Vg (r , t) is the superfluid 
velocity and v„(r, t) the velocity of the normal fluid. For 
^Bcs Rtf, the extent of the vortex core is small com- 
pared to the size of the cloud, and the main effect of 
the vortex on the collective mode spectrum is the pres- 
ence of the vortex velocity field v^, (r) ~ Kh/2maP- 
In these considerations we keep the winding number 
K explicitly different from one for generality. Writing 
n{r,t) = no(r) -|- Sn{r,t) and Vs(r,i) = v„(r) + u(r,t), 
where no(r) is the equilibrium density profile with the 
vortex alone (which we take to be coincident with the 
Thomas- Fermi one without vortex), and linearizing in 
Sn{r,t) and u{r,t), Eq. (UTt) can be written as 
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(18) 



Here, uj is the frequency of the collective mode, -Po(r) 
is the equilibrium pressure profile, and m is the mag- 
netic quantum number of the mode. The velocity field 
associated with the mode has been written as u(r, t) — 
V<I>(r, i) with<I>(r, i) ^ ^{r,9)exp[i{m4>-ujt)]. The term 
Kfim/2maP^ in Eq. ( p^ ) comes from the presence of the 
vortex velocity field v„. Without this term, Eq. (|8|) has 
been solved for a spherical symmetric trap by writing 
$nim( r) = ^ni{r)Yiyn{e,(j)) yie lding the spectrum Unia = 
2wTV'(n2 + 2n + ln + 3Z/4)/3 with n = 0,1,2,... [|6). 
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From Eq. (y_8[), we see that the frequency shift of a given 
mode induced by the vortex can be calculated perturba- 
tively as 
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Here 



*«im|/(r)|*«im) denotes the spatial average 
J^ ^ w{r)dr J dQ^^^ij^{r)f{r) with the weight function 
w{r) = r^(l — r^/i?^^)^/^. This anomalous weight has to 
be introduced in place of simply because the operator 
in Eq. ( p^ without the perturbation is not Hermitian. 

As pointed out in Ref . [ p3[ , the perturbative procedure 
works for |m| > 2; for |r7rp< 2 it predicts an unphysical 
p —^ divergence in the density fluctuation of the mode. 
With no vortex present, the lowest mode for a given angu- 
lar momentum is the surface mode ^n=oim{f) oc with 
frequency VJujt- Recalling that p = rsinO and using the 
fact that (47r)-i/ dl]|yi™(17)f/ sin^ 6* = {21 + 1)/2|to|, 
the matrix elements in Eq. ( [l9| ) can be calculated ana- 
lytically for these surface modes and we obtain for the 
frequency shift: 
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with I TO I > 2. As expected, the vortex splits the 21 + 1 
degenerate modes depending on the direction of the pro- 
jection of their angular momentum on the z-axis. Not 
all the modes are split however since the splitting is 
independent of Jto| in analogy with the equivalent re- 
sult for bosons |23| . Particularly important is the result 
(w| ±2~'^2o)/'^20 — ±V2K/{6Nay^^ for the quadrupolar 
mode I = 2, m = ±2, since this mode is easily excited 
in trapped gases and has been already employed for a 
precise determination of the critical frequency for vortex 
nucleation in Bose gases. 

The same result can be obtained following the sum rule 
approach of ref. j2|] . From that the splitting is found to 
be given by 0^2,2 — ^2,-2 = 2 < Iz > /{ma < >), 
with < Iz > expectation value of the angular momentum 
along the z-axis per atom {k,Ti/2 in the case of a vortex), 
and < p^ > expectation value oix^ +y'^ (equal to Rtf/^ 
for an isotropic cloud with a Thomas-Fermi density pro- 
file) . From the latter result one can immediately see that 
the splitting of the modes of a Fermi superfluid is in gen- 
eral smaller compared with the BEC case, the reason 
being that given the same number of atoms the radius 
of a fcrmionic cloud is usually larger due to the Pauli 
repulsion and thus the expectation value of is also 
correspondingly larger, and the splitting reduced. For 
2Ncr = 10^ particles trapped, the to = ±2 quadrupole 
modes are split by ^ 1%. Although this is a rather small 
shift, it should be measurable assuming the same high 
spectroscopic precision demonstrated for EEC's can be 
obtained for trapped Fermi gases [2^. 



V. CONCLUSION 

In this paper we considered various aspects of the vor- 
tex state of a dilute superfluid Fermi gas at T = 0. For a 
trapped system, we found that a large number of particles 
and a not too small scattering length yields ^bcs Rtf 
and the vortex is well confined within the gas. We then 
used a simple model to calculate the energy of a vortex 
in a uniform medium. Subsequently, using the fact that 
the structure of the vortex near the rotation axis is es- 
sentially unaffected by the trapping potential we derived 
an expression for the vortex energy in a trap, and we 
employed this energy expression to calculate the thermo- 
dynamic critical rotation frequency for the formation of 
a vortex. Finally, we suggested a way of observing the 
presence of the vortex by calculating perturbatively its 
influence on the collective mode spectrum of the gas. In 
the Appendix we report an alternative, less naive, de- 
scription of the vortex in a uniform medium and find a 
slightly different value for its energy compared the one 
obtained in Sec. |^ 
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FIG. 1. Critical number of atoms per spin species for which 
iscs / Rtf = 1 in an isotropic trap. Well above the line the 
local density approximation, and thus Eq. (Q), applies, below 
the line the system is intrinsically finite sized. 



APPENDIX A: A GINZBURG-LANDAU 
DESCRIPTION OF THE VORTEX CORE 

We here present a Ginzburg-Landau description of the 
vortex core in a uniform gas, and the consequent result 
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for the total energy of a vortex in a cylindrical bucket of 
radius Re- As is well-known, Ginzburg-Landau theory 
is only valid for temperatures such that \T — Tc\/Tc <C 1 
but the following calculation can be used for a qualitative 
estimate at T = [||. 

The extension of the Ginzburg-Landau theory to zero 
temperature for a uniform system can be done by impos- 
ing to the free energy [fsl 
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to be equal to the condensation energy density, which 
in a uniform system is given by Cconrf = — 3A§rio-/4eF 
[ p2[ . -0 is here the well known Gin zbu rg-Landau order 
parameter. Upon minimization of (Al) with respect to 
ip* we obtain the Ginzburg-Landau equation 
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For a uniform system the solution is IV'oP = —A/B 
and from the fact that the Ginzburg-Landau free en- 
ergy then coincides with the condensation one fcL = 
-SAlncr/Aep, we obtain A = -3Al/2eF and B = 
-Ajua. 

We now calculate the structure and the energy of the 
vortex. A vortex along the z-axis is described by writing 
the order parameter in cylindrical coordinates as ■i/'f r) — 
f{p)e^'^'^. Replacing this expression into Eq. ( |A2| ) we 
find 
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where we introduced the dimensionless quantities x — 
//IV'ol and x = pj^cL- We used the fact that / does 
not vary along the z direction if the system is uniform, 
and we defined the Ginzburg-Landau coherence length 
^GL = jAmaA, which imphes ^gl = ?iw_F/2\/3Ao = 
0.907^BC75. This equation has exactly the same form as 
the Gross-Pitaevskii equation for a vortex in a uniform 
boson cloud 28|. It can be solved numerically and the 
results for the lowest k (k = 1, 2, 3) was first obtained by 
Ginzburg and Pitaevskii [p^ . A very good approximate 
solution for k, — \ can be obtained by a variational calcu- 
lation yi eldi ng x — a;/(2-|-x^)^/^ Using this solution 
in Eq. (Al), one finds that the energy cost associated 
with the vortex core is given by 



n- Urr In 



2m„ 



\ V2 £.GL 



(A4) 



This result is now identical with that for a Bose- 
Einstein condensate, with the mass of a single bosonic 
atom replaced by that of a Cooper pair {2ma)^ and 



the boson coherence length replaced by the Ginzburg- 
Landau one. Using S^ql = 0.907^505, we obtain = 
TTh^na\ii{D Rc/£,BCs)/'2-ma with D = 1.65. 



[9: 

[lo; 

[11 
[12; 

[is: 

[14 

[is: 
[16: 
[17: 
[is: 

[19 
[20 

[21 
[22: 

[23: 

[24: 

[25: 



[1] B. DeMarco S. B. Papp, and D. S. Jin, cond 



mat/0101445t B. DeMarco and D. S. Jin, Science 285, 
1703 (1999). 

A. G. Truscott et al, Science 291, 2 570 (2001). 



F. Schreck et al, cond- mat/001 1291 



K. M. O'Hara et al, Phys. Rev. Lett. 85, 2092 (2000). 
H. T. C. Stoof, M. Houbiers, C. A. Sackett, and R. G. 
Hulet, Phys. Rev. Lett. 76, 10 (1996). 
F. Weig and W. Zwerger, Europhys. Lett. 49, 282 (2000); 
W. Zhang, C.A. Sackett, and R.G. Hulet, Phys. Rev. A 
60, 504 (1999); J. Ruostekoski, Phys. Rev. A 60, R1775 

(1999) ; P. Torma and P. Zoller, Phys. Rev. Lett. 85, 487 

(2000) . 

See e.g. A. L. Fetter and A. A. Svidzinsky, J. Phys. Cond. 
Mat. 13, R135 (2001). 

M. Rodriguez , G.-S. Paraoanu, and P. Torma, cond 



mat/0104491 



R. I. Epstein and G. Baym, Ap. J. 328, 680 (1988). 

L. P. Gorkov and T. K. Melik-Barkhudarov, Sov. Phys. 

JETP 13, 1018 (1961). 

H. Heiselberg, C. J. Pethick, H. Smith, and L. Viverit, 
Phys. Rev. Lett. 85, 2418 (2000). 

P. G. de Gennes, Superconductivity of Metals and Alloys 
(Benjamin, New York, 1966). 

A. L. Fetter and J. D. Walecka, Quantum Theory of 
Many-Particle Systems (McGraw-Hill, New York, 1971). 
G. Bruun, Y. Castin, R. Dum, and K. Burnett, Eur. 
Phys. J. D 7, 433 (1999). 

R. Combescot, Phys. Rev. Lett. 83, 3766 (1999); H. 
Heiselberg, Phys. Rev. A 63, 043606 (2001). 
P. Nozieres and D. Pines, The Theory of Quantum Liq- 
uids (Benjamin, New York, 1966), vol H. 
V. L. Ginzburg and L. P. Pitaevskii, Sov. Phys. JETP 
34, 858 (1958). 

E. Lundh, C. J. Pethick, and H. Smith, Phys. Rev. A 55, 
2126 (1999). 

F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999). 

P. Torma and P. Zoller, Phys. Rev. Lett. 85, 487 (2000); 

G. M. Bruun, P. T orma, M. Rodriguez, and P. Zoller, 
cond-mat/00 11333 . 

F. ZambeUi and S. Stringari, Phys. Rev. Lett. 81, 1754 
(1998). 

S. Sinha, Phys. Rev. A 55, 4325 (1997); R. J. Dodd, K. 
Burnett, M. Edwards, and C. W. Clark, Phys. Rev. A 
56, 587 (1997). 

A. A. Svidzinsky and A. Fetter, Phys. Rev. A 58, 3168 
(1998). 

F. Chevy, K. W. Madison, and J. Dahbard, Phys. Rev. 
Lett. 85, 2223 (2000). 

L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perg- 



6 



amon Press, Oxford, 1987). 
[26] G. M. Bruun and C. W. Clark, Phys. Rev. Lett. 83, 5418 

(1999); M. A. Baranov and D. S. Petrov, Phys. Rev. A 

62, 041601(R) (2000). 
[27] D. S. Jin et al, Phys. Rev. Lett. 77, 420 (1996); D. M. 

Stamper-Kurn et al, Phys. Rev. Lett. 81, 500 (1998). 
[28] C. J. Pethick and H. Smith, Bose-Einstein Condensation 

in Dilute Gases (to be published by Cambridge University 

Press). 

[29] A. L. Fetter, In Lectures in Theoretical Physics, eds. K. 
T. Mahanthappa and W. E. Britten (Gordon and Breach, 
New York, 1969), Vol. XIB p. 351. 



7 



